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The correlated basis function theory is applied to the study of medium-heavy doubly closed 
shell nuclei with different wave functions for protons and neutrons and in the jj coupling scheme. 
State dependent correlations including tensor correlations are used. Realistic two-body interactions 
of Argonne and Urbana type, together with three-body interactions have been used to calculate 
ground state energies and density distributions of the 12 C , ie O . 40 Ca , 48 Ca and 208 Pb nuclei. 
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PACS numbers: 21.60.-n,25.30.Dh 

I. INTRODUCTION 

In the last decade, the validity of the non relativistic many-body theory in the description of nuclear systems, has 
been well established. The idea is to describe the nucleus as a set of mutually interacting nucleons with a hamiltonian 
• of the type 

C2> — h 2 

\o~. # = E ^- v * + E «« + E "«* (!) 

Z~_~j • i ' i<j i<j<k 

where the two- and three-body interactions, and Vijk respectively, are fixed to reproduce the properties of the two- 
and three-body systems. 

Quantum Monte Carlo calculations solve the many-body Schrodinger equation without approximations 0, Q and 
have been applied with success to describe nuclei up to twelve nucleons Q. The only limitation of this approach is 
related to the available computing power, which, at present, does not allow its straightforward application to medium 
and heavy nuclei. The good results obtained by the Quantum Monte Carlo calculations have renewed the interest 
in nuclear structure calculations with realistic, or microscopic, interactions. For this reason, many-body techniques, 
which are not approximation free, have been formulated, or reconsidered, and applied to the study of medium and 
heavy nuclei, or more in general, to nuclear systems with A>4. 

The oldest, and most commonly used, technique is the Brueckner Hartree-Fock approach recently applied to cal- 
culate the properties of 16 O 0,0. In the last years, no-core shell model calculations 0,0 have used microscopic 
interactions to calculate some property of nuclei up to 12 C . The coupled cluster method has been used to evaluate 
16 O properties HHE3- 

More recently a technique called Unitary Correlation Operator Method has been developed 
to produce effective interactions used in mean-field and shell model calculations, and it has been applied with success 
to a variety of nuclei up to 208 Pb [Tllll2l|. 

In the early nineties, we started a project aimed to apply to finite nuclei the Correlated Basis Function (CBF) 
theory, originally developed for infinite systems |l3|. The idea was to use the Fermi Hypernetted Chain (FHNC) 
resummation technique, well tested in nuclear matter, to describe finite systems. We first tested the validity of our 
approach in model nuclei with degenerate proton and neutron states, and in Is coupling scheme |l ll Il5|. Then, by 
using central interactions and scalar correlation functions, we applied the theory to describe nuclei with different 
single particle bases for protons and neutrons, with different number of protons and neutrons and with the more 
realistic jj coupling scheme |16). More recently, we have treated two-body interactions containing tensor and spin- 
orbit terms and this implies the use of state dependent correlations |l7^ . This was not a straightforward extension 
of the formalism because the non commutativity of the different operators requires the use of the Single Operator 
Chain (SOC) approximation in the FHNC resummation. The FHNC/SOC computational scheme has been extended 
to include spin-orbit and three-body interactions of Urbana and Argonne type |l8j . In these last calculations protons 
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and neutrons single particle wave functions were distinguished only by the third component of the isospin, and they 
were treated in the Is coupling scheme. This choice allowed us to treat only systems saturated in spin and isospin. 

This paper represents a step forward in our project. The FHNC/SOC formalism has been extended to treat nuclei 
with different proton and neutrons wave functions, and in a jj coupling scheme in spherical basis. With this single 
particle basis we obtained results for the 12 C , 16 , 40 Ca , 48 Ca and 208 Pb nuclei. 

The aim of this paper is to present the results we have obtained by using the two-body Argonne v' s 0, El and 
the Urbana V14 |2Ct l2lj interactions, where we have considered all the channels up to the spin-orbit ones. In our 
calculations we have implemented these two-body interactions with the Urbana IX [j| and the Urbana VII [2^ 
three-body interactions respectively. 

In Sect. |H]we briefly discuss the basic ideas of the formalism pointing out the main differences with respect to 
what has been presented in and 0. The detailed presentation of the heavy formalism, which is given in j2^|, is 
beyond the scope of this article. 



II. THE FORMALISM 



Our approach is based on the variational principle, i.e. we search for the minimum of the energy functional 

Em = < ™* > 



(2) 



We used the hamiltonian of eq. m our calculations, and we express the two-body interaction in the form 

v ij = v p( r ij)°^ ( 3 ) 



where we have defined the operators O p as: 



Ofr 1 ' 8 = [1, cr, ■ cr^Sv, (L • S) -] ® [l, Ti • Tj] 



and where we have indicated with Sij the tensor operator 



Sij — 3- 



- cr, ; 



(4) 



(5) 



In eq. (J3J) we have limited the sum to the first 8 channels which are those used in our calculations. In reality modern 
two-body interactions have been parametrized by using up to 18 different channels fl9j |. 

The three-body interaction |2Ll22j| is composed by two parts. A first one describes the process where two pions are 
exchanged with the intermediate excitation of a A j24[. This term describes the long-range part of the three-body 
interaction and, usually, gives an attractive contribution. The second part of the three-body interaction is a shorter 
range, repulsive, spin and isospin independent term. This term is supposed to simulate dispersive effects obtained 
when the A degrees of freedom are integrated out. 

The search for the energy minimum (J2J is done in the subspace of the many-body wave functions $ which can be 
expressed as: 



*(1,2,...,A) = F(1,2,...,^)$(1,2,...., J 4) 



(6) 



where F(l,2, A) is a many-body correlation function and $(1, 2, A) is a Slater determinant composed by a 
set of single particle wave functions 4> a {i)- In our calculations we describe the many-body correlation function as a 
symmetrized product of two-body correlation functions: 



F(1,2,...,A)=S 



(7) 



and we use two-body correlation functions Fij which have an operator dependence analogous to that of the nucleon- 
nucleon interaction ©, but without the spin-orbit terms: 



p=l,6 



(8) 
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A. The cluster expansion 



The evaluation of the multidimensional integrals necessary to calculate the energy functional (J2J) is done by per- 
forming a cluster expansion and by applying the FHNC resummation techniques. Our work is based on the formalism 
developed in [3 El HHH. 

The two basic quantities used in the cluster expansion of eq. (01 are the one- and the two-body density matrices 
respectively defined as: 

Pxixux'-t) = A f dx 2 dx 3 ...dx A y*(x 1 ,x 2 ,...,x A )y(x' 1 ,X2,...,x A )/ < *|* > (9) 

P2(XI,X 2 ) = A{A-1) J dx 3 ...dxA^*(x 1 ,X2,X 3 ,...,XA)^(xi,X2,X 3 ,...,X A )/ < > (10) 

where we have indicate with Xi the coordinate fully characterizing the particle i, i.e. position r^, 3-rd components of 
the spin, s, and of the isospin, t. The evaluation of the expectation values of one- and two-body operators can be 
obtained by making the operators acting on the densities © and (|10|) and integrating on the free variables. 
In the uncorrelated case, F(l, 2, A) = 1, the expressions © and (fTTTft become: 

Plixx,^) = ^(xr^K) (11) 

i 

p 2 {x u x 2 ) = p°(xx)p°(x 2 ) - p 1 {x 1 ,x 2 )Pi(x2,x 1 ) (12) 

where we have indicated with </>j the single particle wave functions. In the above equations, the sum is done on all 
the occupied states, and we have used the notation p°(xi) = p1{x\, x\) to indicate the diagonal part of the one-body 
density function. 

We construct the single particle wave functions by using a spherical mean field potential containing a spin-orbit 
term. Protons and neutrons wave functions are generated by different potentials. The single particle wave functions 
are classified in terms of the angular momentum j and of the third component of the isospin t. They can be written 
as: 

^nijm^i) = Rt mj(n)J2 < l ^\ s \] m > ^(fii)x.(*)Xt(») (13) 

where we have indicated with Yi^(£li) the spherical harmonics, with Xs,t{i) the Pauli spinor/isospinor and with 
< lp^s\jm > the Clebsch-Gordan coefficient. 

In the calculations of Refs. 0, 0, 0], the radial part of the single particle wave function was independent 
on the isospin. This, joined to the fact that we studied doubly closed shell nuclei on Is coupling with the same 
number of protons and neutrons, allowed us to separate the spatial part from the spin and isospin dependent parts in 
p°{xi, x[). Because we restricted our investigation to finite systems saturated in spin and isospin, we could use the 
trace techniques developed in [25| to deal with state dependent correlations. 

In the present case, as pointed out in |l6|, the choice of the single particle basis l|13|l requires to consider the spin 
dependence of the radial part of the one-body density matrix: 

plixux'J = J2 pf\r l ,r' 1 )xi(l)xAl')xt(l)Xt(l') (14) 

s,s' ,t 

therefore the factorization between spatial and spin-isospin dependence of the one-body density matrix is lost. 

We consider two radial terms of the one-body density matrix. A first one, Pq, related to the parallel spin terms, 
and the other one, p^-, to the antiparallel spin terms. We can rewrite eq. (fH|) as: 



x[) xt (i)xt(i') I PhinA) ^ xl(i)x.(i') + PoMA) ^(-ir i/2 xJ(i)x- s (i') ) (15) 

t \ s s 

The explicit expressions of p l Q and pl,^ in terms of single particle wave functions, are given in 16]. 
The study of Ref. 

has shown that the main contribution to the one-body density matrix is given by p^, which 
is the only term surviving in spin saturated systems having degenerate spin-orbit partners. The similarity of this 
term to the one appearing in the Is coupling scheme, allows us to use some simplifications in the calculations of the 
correlated distributions I|9I10|I following the procedure adopted in [2^| . Even though the contribution of the terms 



(l)Xt(l') (Wi, ri) £ xt(l)Xs(l') + pl 3 {viA) E^ 1 )^ 17 ^ (l)x-(l') 

\ s s / 
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are much smaller than those of p* Q , they cannot be neglected, as it is indicated by the sum rule results. For this reason, 
we had to calculate new kinds of spin traces. 

In the calculation of the one-body density matrix and of the two-body density function and (|lt)(l we consider the 
scalar terms of the correlations at all the orders of the cluster expansion. This part of the calculation follows step by 
step that presented in |l6j . The two-body densities are written in terms of the scalar correlation fi(r), of the nodal, 
or chain, diagrams N xy (ri, r 2 ), of the non nodal (composite) diagrams X xy (ri,T2) and of the elementary, or bridge, 
diagrams E xy (v\, V2). Here the two sub- indexes x and y identify the exchange character of the external points 1 and 
2, as classified in [lij. 

The situation becomes more complicated when state dependent correlations are introduced, essentially because the 
various terms of the correlations do not commute with the hamiltonian and among themselves. At present, a complete 
FHNC treatment for the full, state dependent cluster expansion of the two-body densities is not available. In our 
calculations we adopted the Single Operator Chain (SOC) approximation |25| consisting in summing those operator 
chains containing scalar dressing at all orders, but only a single operator p > 1 element. 

The different kind of diagrams: nodal, non-nodal and elementary ones, are characterized by an index p denoting 
the operator (jJJ associated to the correlation function. In order to treat nuclei with different number of protons and 
neutrons (N ^ Z) we separate the isospin dependence of the correlation operators from that of the other operators: 

F ^l) = E E hk-i + i(r l3 )0^- 1+l = J>, • Tj ) 1 ]T hk-x^rn)^ (16) 

fc=l (=0 1=0 k=l 

with pk -1 ' 2 ' 3 = 1, cri ■ o-j, Sij. 

With respect to the treatment presented in 0] , now we have to consider that the various diagrams depend on the 
third isospin component of the nucleons located in the two external points. As a consequence, the number of cluster 
diagrams to be calculated is four that of the previous cases. In reality, for symmetry reasons, we had to calculate 
separately only three different types of diagrams. In addition, we should take into account also the contribution of 
the uncorrelated two-body densities with antiparallel spins, the terms related to p\j in eq. I|15|) . The contributions 
of these two-body densities are not zero only for the scalar, p = 1, and isospin, p — 2, operators. As we have already 
mentioned in the introduction, a detailed presentation of all the expressions used in our calculations can be found in 

m 



B. The energy expectation value 

The calculation of the energy expectation value has been done along the lines of Refs. and 01- We used the 
Jackson- Feenberg identity to calculate the kinetic energy, since this expression allows us to eliminate terms of the type 
(VjF)(Vj$) involving three-body operators (2(|. Following the notation of Ref. [14|. we express the kinetic energy 
by separating those terms where the V operator acts on the correlation function F only, from those where it acts on 
the uncorrelated many-body wave function $. 

* = T JF =T <t> +T F , (17) 



where we have defined \l< 



and 



h 2 

Ts = -A — ($*T 2 V 2 $ - (Vi$*) F 2 (Vi$)) (18) 
4m 

h 2 

Tp = —A — ($* \FVlF- VxF- Vi.Fl (19) 
4m 

The operator structure of the Tf terms is similar to that obtained by inserting the two-body interaction term of the 
hamiltonian. For this reason, we found convenient to calculate the joint contribution of the Tp + V2 = W operator, 
called interaction energy in 17]. In this calculation we have considered all the interaction terms up to p = 6. The 
expectation value of W can be expressed in terms of 



H p jf r (r 12 ) ee W^ r (r 12 ) = T^'> 12 ) + V^ r {r 12 ) 

' ■ (~ 2Vi Sql {/ P ( r i2)V 2 r (ri 2 ) - V/^(r 12 ) • Vf> 12 )} + P(r 12 )v* (r 12 )f (r 13 )) 



(20) 
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where the p, g, r labels refer to the type of operator and can assume values from 1 to 6. Following what has been done 
in P"t||. we found convenient to split W in four parts identified by the different topology of the various links, 



W = Wo + W s + W c + W c . 



(21) 



The Wo term is the sum of the diagrams having only scalar chains between the interacting points, connected by Hjp- 
We indicate with W s the sum of the diagrams with a Single Operator Ring (SOR) touching the interacting points and 
scalar chains. The W c term contains diagrams with a Single Operator Chain (SOC) between the interacting points 
and W cs contains both SOR and SOC between the interacting points. The expressions we obtain for the contributions 
of the various W terms are more involved than those presented in |17| . This because we have to sum on all the isospin 
indexes and, in addition, we must consider also the contribution of the exchange links with different spin projections. 
The detailed expressions are given in |2^|. In reality, we calculate explicitly the Wq, W s and W c contributions, while 
for the W cs we used the approximation W cs ~ WcfW^/Wo] whose validity is discussed in [171 ]. 

As an example of the differences between the present calculations and those of 01 j we compare the expression of 
Wq in Is coupling scheme and N = Z that is 



Wo = ~ 



\ J dv 1 dv 2 lP J f r {n 2 ) 



K pqr -A r >j,dir(l, 2) + A n K n r m Ki rm A m P2 , exc (l, 2) 



(22) 



with the present expression of Wo given by 



1 



W = - I dr 1 dr 2 H JF 



2fci — 1+Zi,2fc 2 — l+/ 2 ,2fc 3 — 1+Z 3 



(ri 2 ) 



1 A P2,diA L ' Z )Xl 1+ l 2+ l 3 + 



( a/3 . a/3 \ a/c 4 rfe 4 fe 1 fe 5 rk 2 k 3 k 6 ( T k 5 k e l a/3 /-, n\ . rk 5 k e 2 afi r>\\ 



(23) 



In both equations, a sum on repeated indexes is understood, with p, q,r,m,n = 1,...,6 in llL'l'l) and a,@ — p,n, 
k{ = 1, 2, 3 and ij = 0, 1 in (|23() . We have indicated with p 2 ,dir the direct part of the two-body density function which 
includes all the diagrams where the two external particles, 1 and 2, do not belong to the same statistical loop. Eq. 
has only one exchange part of the two-body density function , p2.exe- On the contrary eq. (|23|l contains two 



terms depending weather the spin of the external particles are parallel, p 2 f exc , or antiparallel, p 2 P exc: 

While the two-body density functions are isospin independent in eq. Q22[l. in the present calculations we have to 
distinguish the cases when the external particles are protons or neutrons. This isospin dependence obliges us to isolate 
the isospin traces from those related to the other operators. These traces can be written in terms of 



a/3 



y a/3 

A 1 1 



X* a (l)x}V) (T 1 -r 2 ) n Xa(l)x/3(2) 



(24) 



In eq. i|22|) the isospin trace is included in the K matrices, defined in [2f|. The I matrices of eq. EUare the spin traces 
of the products of operators, therefore, they are included in K. As we have already mentioned, the presence of a new 
statistical link in the exchange part of eq. I|23|l requires the evaluation of spin traces absent in the Is coupling case. 

The differences pointed out in the previous discussion are present in all the terms of the kinetic energy. To simplify 
the calculations, we have calculated the contribution of p 2 ^ exc j only in the Wo and W s parts of the interaction energy. 

In these two terms the contribution of p 2 exc j is about two orders of magnitude smaller than that of p 2 exc - We have 
estimated that this relation between the two exchange terms is not going to be modified in W c . 

The calculation of the kinetic energy contribution requires the evaluation of the term of eq. I|18|) . Also in this 
case, in analogy to what we have done in |l6lll7| . we separate the term in three parts: 



Tj, — 



,{2) 



T 



(3) 



(25) 



where the term n = 1 corresponds to the contributions where the kinetic energy operator acts on a nucleon not 
involved in any exchange, the term n = 2 when it belongs to a two-body exchange loop, and n = 3 to a many-body 
exchange loop. 

The expression obtained for is the same as that given in [l7| . The other terms are more involved because of 
their dependence on the isospin indexes, the presence of new isospin traces and of the exchange links with the p° 
uncorrelated densities |23j . 

The evaluation of the interaction energy l|2l)|l . allows us to calculate all the two-body interaction channels up to 
p = 6. The contribution of the remaining terms of the two-body interaction, the p = 7, 8 spin-orbit channels, has been 
calculated following the procedure and the approximations discussed in ■ We calculate only the Wq term with the 
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fe correlation function. Again, the expressions we have obtained are more involved than those given in [18( since we 
have isospin dependent terms, and we have considered the contribution of the leading p° in the exchange terms. 

The above discussion can be extended to describe our treatment of the three-body interaction. We have used the 
nuclear matter results of Ref. [27| to select the relevant diagrams to be evaluated. We followed the lines of calculations 
presented in and we took care of the new aspects caused by the jj coupling and N ^ Z |2^| . 

III. RESULTS 

The FHNC/SOC formalism presented above, has been applied to calculate binding energies and matter distributions 
of 12 C , 16 , 40 Ca , 48 Ca , and 208 Pb doubly magic nuclei. A complete variational calculation would imply the search 
for the energy minimum by modifying both the set of single particle wave functions forming and the correlations 
function F. These calculations are numerically very involved since the number of the free parameters to be modified 
is rather large, about ten. To reduce the computational effort in our calculations we kept fixed the single particle wave 
functions and we made variations on the correlation function only. At first sight this can appear a severe limitation, 
however we have shown in jl8l | that large changes of the mean field wave functions produce small variations of the 
value of the total energy of the system. For this reason, by extrapolating the findings of 01 we are confident that 
our results are only few percent above the minimum obtained with a full minimization. 

In our calculations we adopted the same set of single particle wave functions used in |Tfij . They are generated by 
a Woods-Saxon central potential containing a spin-orbit term. Different potentials have been used for protons and 
neutrons. The parameters of these potentials, given in [16j, have been taken from the literature where they have been 
fixed to reproduce, at best, the single particle energies close to the Fermi energy, and the charge density distributions. 

The other fixed input of our calculations are the nucleon-nucleon interactions. In this article we discuss mainly the 
results obtained with the two-body Argonne v' B interaction 0,0] implemented with the Urbana IX {UIX) three-body 
force . As a comparison, we discuss also results obtained with the Urbana two-body interaction [23, El but 
considering o nly the first eight channels. Together with the V14 two-body interaction interaction we used the Urbana 
VII (UVII) H3 three-body interaction. 

The hamiltonian and the single particle basis are fixed input of our calculations and the variations have been done 
on the correlation function F. The two-body correlation functions f p (r) are the functions we have modified in order 
to find the minimum of the energy functional. We fixed these functions by using the technique we have defined as 
Euler procedure in our previous works 0,^3- We solve the FHNC/SOC equations up to the second order in the 
cluster expansion by forcing the asymptotic behaviour of the two-body correlation functions f p (r). More specifically 
we imposed that, after certain asymptotic distances, that we call healing distances, the scalar two-body correlations 
(p = 1) must be equal to one, while the other correlation functions must be zero. The values of the healing distances 
are our variational parameters. 

In the actual calculations, in order to reduce the number of variational parameters, we used only two healing 
distances. A first one, d c , for all the central channels, p < 4, and a second one, dt, for the two tensor channels 
p = 5, 6. Nuclear matter [2f| and finite nuclei [T^. fl8| calculations indicate that tensor and central correlations heals 
at different distances. Minima are found for tensor healing distances larger than the central ones. 

In Tab. [IJ we present the values of the healing distances giving the energy minima for each nucleus, and for both 
interactions v' s with UIX and V14 with UVII. It is remarkable the similarity of these values for all the nuclei but 
the 12 C nucleus. The values of d c for the two interactions are the same, with the usual exception of 12 C , while the 
dt values for U14 plus UVII are slightly larger than those obtained with v' 8 plus UIX. 

The two-body correlation functions obtained for the v' s plus UIX interaction are shown in Fig. ^for each channel. 
In the figure, we also show with the diamonds the nuclear matter correlations. These last correlations are obtained 
with the Euler procedure by using the same healing distances used for finite nuclei but with constant density value of 
the nuclear matter the saturation point. Analogous results have been obtained for the U14 plus UVII interaction. 

We would like to point out two features of these results. The first one, is that the various correlations are rather 
similar for all the nuclei considered. Only the 12 C results are out of this general trend. This fact indicates that 
the correlation functions are more sensitive to the characteristics of the nucleon-nucleon interaction than to the shell 
effects. The second point we want to mention is the difference between the healing distances of the central and tensor 
channels. In the two-body interaction, the tensor channels are active at larger distances with respect to the central 
channels, and this produces the consequences on the correlation functions we have just pointed out. It is interesting 
to notice that the healing for the tensor channels appears at distances larger than the experimental charge radius of 
the deuteron. 

In this discussion the 12 C nucleus remains an anomaly. The fi correlation is rather similar to those of the other 
nuclei, all the other two-body correlations are remarkably different, especially f% and f± related to the spin and 
spin-orbit dependent terms. 
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Also the nuclear matter results are not on the band formed by the correlations found for the 16 O - 208 Pb nuclei. 
The differences are not so large as those of the 12 C , but still rather evident, ft is remarkable the similarity in the /i 
channel and in the fs channel, while all the other channels show noticeable differences. 

Our calculations are organized in such a way that we first solve the FHNC equation for the scalar correlation only, 
and in a second step we solve the FHNC/SOC equations with the other state dependent correlation terms. The 
validity of our solution techniques, from both the theoretical and numerical points of view, can be appreciated by 
observing the convergence in the density sum rules shown in Tab. [HI In this table we give the normalization of 
the one- and two-body densities © and 1|1U|) . The Si and S2 sum rules have been defined in such a way that all 
of them are normalized to unity. The results of the table show that the Jastrow-FHNC calculations, those with fx 
only, produce densities deviating by few parts on a thousand from the correct normalization. The inclusion of state 
dependent terms in the correlations, the FHNC/SOC calculations, is responsible for an increase of these deviations 
up to a maximum value of 5.5 %. In many-body systems saturated in spin and isospin one can also test spin and 
isospin density sum rules. In the systems we are investigating these sum rules are not any more valid. 

The main result of our paper is summarized in Tab. IIIII where we give the values of the binding energies per nucleon 
for all the five nuclei considered, and we compare them with the experimental binding energies [28j . In the table we 
present the various terms contributing to the total energy: the kinetic energy T, the two-body interaction, where the 
the contribution of the first six channels V 2 6 _ body an d that of the spin-orbit interaction Vls are separately given, the 
Coulomb interaction Vcoui and the three-body force V-i-body In the kinetic energy term the contribution of the center 
of mass motion, calculated as discussed in has already been subtracted. 

The various terms show some saturation properties. For example the values of kinetic energies per nucleon T 
increase up to 40 Ca and then they remain almost stable around a value of about 40 MeV. An analogous behaviour is 
shown by the V^ 6 _ bod terms. The contribution per nucleon increases with increasing number of nucleons up to 40 Ca 
. After that it is evident that saturation appears in heavier systems. 

We have already mentioned the fact that the spin-orbit terms are not treated consistently in the FHNC/SOC 
computational scheme, but they are evaluated by using some approximation. In any case, in all the nuclei considered, 
their contributions are of the order of few percent with respect to the V' 2 6 _ bodj( contributions. As a further test, we 
have done calculations in 16 and 40 Ca switching off the spin-orbit terms in the mean field potential. In this case the 
spin-orbit partner single particle wave functions are identical. The differences in the total spin-orbit contributions, 
with respect to the values given in Tab. IIIII are within the numerical uncertainty. 

The contribution of the Coulomb term Vcoui is evaluated within the complete FHNC/SOC computational scheme. 
As expected, the behaviour with increasing size of the nucleus does not show saturation because of the long range 
nature of the interaction. The Coulomb terms behave as expected, their contributions increase with increasing number 
of protons. The apparent inversion of this trend from 40 Ca to 48 Ca is due to the representation in terms of energy 
per nucleon, which in this case is misleading, since the proton number is the same for the two nuclei. In this case it 
is better to compare the total values of the Coulomb energies, 78.40 MeV for 40 Ca and 75.36 for 48 Ca . The 3.8% 
difference between these to values is due to the different structure of the two nuclei. The inclusion of the Coulomb 
repulsion reduces the nuclear binding energies. To simplify the reading of the table we give in the raw labeled T + V(2) 
the intermediate energy value obtained with the two-body force only, Coulomb interaction included. The total binding 
energies without the three-body forces are all above the experimental values. 

In addition to that, there is the contribution of the three-body force. As discussed in [l^, the two terms composing 
this interactions provide contributions of different sign, the Fujita-Miyazawa term is attractive, while the other term 
is repulsive. In our calculations the total contribution of the three-body forces is always globally repulsive. 

The comparison with the experimental energies indicates a general underbinding of about 3-4 MeV per nucleon. In 
the general behaviours we have discussed so far, the anomaly of the 12 C nucleus emerges. In all our calculations the 
attractive contribution of the V^ 6 _ body is relatively smaller than for the other nuclei. As a consequence, this nucleus 
results barely bound in our calculations. Some crucial physics ingredient relevant in 12 C while negligible for the other 
other nuclei is missing in our approach. Probably this has to do with soft deformations of the 12 C nucleus which we 
are unable to treat. 

The comparison between the two interactions indicates that the vu plus UVII is more attractive than the v' s plus 
UIX. This fact is already present when the two-body interaction only is considered. The situation is enhanced by 
the three-body force, more repulsive for the UIX case than for the UVII one. The contributions of the spin-orbit 
term in the two cases have different sign, they are attractive for v' s and slightly repulsive for vu. The differences in 
the total energies calculated with the two sets of interactions go from a minimum of 5% ( 16 ) to a maximum of 18% 
( 208 Pb). 

Despite the remarkable differences in the energy, the two interactions produce very similar results in the observables 
related to the density distributions. For this reason, henceforth, we shall present results obtained with the v' s plus 
UIX interaction only. 

In Figs. |21 and |3| we show the one-body density distributions for protons and neutrons respectively. The full lines 
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show the uncorrelated distributions, the dotted lines those obtained by using scalar correlations only /i, we call 
Jastrow these calculations, and the dashed lines the results of the complete calculation. The Jastrow results produce 
distributions which are smaller at the center of the nucleus with respect to the mean-field distributions. This effect is 
strongly reduced when all the correlations are included in the calculation. These findings are in agreement with the 
results of ref. .29] where a first-order cluster expansion has been used. It is evident the scarce relevance of short-range 
correlations in reducing the occupation probability of the 3sl/2 proton state in 208 Pb . As pointed out in [^3, the 
effect of long-range correlations seems to be main reduction mechanism. 

In order to investigate the relative effects of the various correlations on the protons and neutrons distributions we 
studied the ratio: 

where, as it has been defined in eq. p°(r) is the uncorrelated, mean-field, density. The Ai ratios are shown in 

Fig. ^ for all the nuclei under investigation. The behaviour of the various lines indicates more similarity between 
results obtained with the same type of correlation than between protons and neutrons. This result is quite remarkable 
because of the large difference between proton and neutron distributions, especially for 48 Ca and 208 Pb . 

A direct comparison with the empirical charge density distributions is not very meaningful since the mean-field 
potential has already been fixed to reproduce, at its best, these distributions. The only statements we can make is 
whether the inclusion of the correlations improves or not the agreement with the data. For this reason we show in 
Fig. the ratios 

A 2 (r) = P(r) "^ (r) (27) 

where the experimental charge distributions have been taken from [3lj . The inclusion of correlations lowers the charge 
density distributions at the center of the nuclei. If scalar correlation only are considered the effect is too large. When 
all the correlations are included the effect is strongly reduced and the agreement with the empirical densities improves. 

The effects of the short-range correlations are better shown on the two-body density matrices l|l(J|l . We define the 
operator dependent two-body density matrices as: 

p * (ri ' r2) ~ 

dx 3 ...dx A ^*(x 1 ,x 2 ,x 3 ,...,x A )V a (l)O p (x 1 ,x 2 )V f3 (2)^(x 1 ,x 2 ,x 3 ,...,x A ) (28) 

where we have indicated with V a (i) the operator selecting protons or neutrons and we have understood the trace 
on the spins and isospins of the particles I and 2. We calculate the values of (|28|l in terms of relative distance 
T\i = |ri — r 2 | between the two considered nucleons 

pf' P (ri2) = J dK 12P f' p '(n.ra) (29) 

where R 12 = (ri + r 2 )/2 is the center of mass coordinate of the pair. In the FHNC/SOC computational scheme, we 
calculated the two-body density matrices as: 

af3,2k 3 -l+l 3( \ _ f2k 1 -l+l 1 {ri 2 )f 2 k 2 -l+l 2 (ri 2 ) J k k k Ak2-,a0 <*P (- -\ (on\ 

P2 (ri,r 2 J- fl{r 12 ) | Xi 1+l2+ i 3 P 2 ^ dir {ri, r 2 J 

if rkikxkn rk 2 k 3 k 5 Ak 5 aP i \ . rfc 4 feife 5 Tk 2 k 3 k(j rfc 5 fc 6 2 a0 / \\ a fc 4 \ a/3 I 

+ { 2 1 A P2,excl r l' r 2j + 1 1 1 P2,excj{ r l> r 2))A X h+h+h+U { 

h=0 ) 

where a sum on the repeated indexes is understood and the symbols of eq. (|23[) have been used. 

Let's consider first the Jastrow case, where only /i correlation functions are present, i.e. 2ki — 1+h = 2k 2 — l+l 2 = 1. 
In this case, apart from a normalization factor, eq. (|28|l gives the probability density of finding two nucleons with 
isospin third-components a and (3 at the distance 7*i 2 . These scalar two-body density functions, or probability densities, 
are shown in Figs. H3 [7\ and 00 for all the nuclei we have investigated. The dotted lines represent the uncorrelated 
joint probability densities of finding the two nucleons as often given in the literature, see for example j^]. This 
joint probabilities are obtained as product of the uncorrelated one-body densities. This definition of the uncorrelated 
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two-body density can be meaningful from the probabilistic point of view, but it is misleading in our framework, since 
it corresponds to use only the first, direct, term of eq. (|12|) . In fermionic systems, the uncorrelated two-body density is 
given by the full expression (|I2[) containing also the exchange term. These complete uncorrelated two-body densities, 
are shown in Figs. El 13 and [S] by the full lines. The effects of the short-range correlations can be deduced by comparing 
these lines with the dashed-dotted lines obtained with the /i scalar correlations only, we call these Jastrow results, 
and with the dashed lines showing the full FHNC/SOC results. 

In all our results, the correlations reduce the values of the two-body densities at short internucleon distances. The 
exchange term of the uncorrelated density, already contributes to this reduction, but the major effect is produced 
by the short-range correlations, and mainly by the scalar correlations. The effects of nucleon-like two-body densities 
are rather similar, as it can be seen comparing Figs. El and El When the pair of nucleons is composed by different 
particles the situation slightly changes. Beside to the strong reductions at small distances, the correlations produce 
enhancements around 2 fm in all the nuclei considered. 

In order to discuss the effects of the correlations on the other operator dependent two-body densities l|28|l . we show 
in Figs. l9l and ITU1 the results for the 208 Pb nucleus only. The results of the other nuclei considered have similar features 
pflj . In these figures, we used the same line conventions of Figs. EI d and [HI When the external particles are nucleons 
of the same type, i.e. a = (3 in eqs. I|28|) - l]3ip. the isospin trace in eq. [2H1 is one - I n this case, the p=l,2,3 two-body 
densities are equal, respectively, to those calculated with p=2,4,6. The same idea can be expressed more formally by 
stating that, in eq. 1|31|) for a = (3 we have P2 a ' 2k ~ 1 (V12) = p^"' 2 * ( r i2)i since Xn° = 1 f° r every value of n. For this 
reason, in Fig. El we present only the p=3 and p=5 two-body densities for both proton-proton and neutron-neutron 
cases. In the proton-neutron case the two-body densities are different in every channel. 

We first notice that tensor two-body densities, p=5, and also p=6 in Fig. [7| are different from zero only in the full 
FHNC/SOC calculations. This occours because in order to get a spin trace different from zero in eq. EH at least two 
tensor operators are needed. In Fig. 0also the p=3 spin density is different from zero only in the FHNC/SOC case. 
The reason is that, in the calculation of this densityfor the uncorrelated and the Jastrow cases, the first term of eq. 
(|I2J) . the direct one, is zero because of the spin trace, while the exchange term is zero because of the different isospin 
third component. 

We observe from Figs. El and that the effects of the correlations are analogous to those already discussed for the 
scalar two-body densities: the values of the densities is strongly reduced at short internucleonic distances. In addition 
to that it is worth to mention the change of sign in the isospin, p=2, densities in the proton-neutron case. The most 
remarkable feature is the range of the various densities. The scalar densities of Figs. and El an d the isospin 
densities, p=2, extend up to relative distances comparable to the dimensions of the nucleus, 15 fm. On the contrary, 
all the other two-bod y d ensities have much smaller ranges, of the order of 3-4 fm. These features are present in all 
the nuclei considered 1231. 



IV. SUMMARY AND CONCLUSIONS 



With this paper we made an important step forward in our project aimed to provide a microscopical description 
of medium-heavy nuclei, by using the CBF theory. With respect to the works of Refs. 0, 0] we have applied for 
the first time the FHNC/SOC computational scheme to nuclei with different number of protons and neutrons and in 
the jj coupling scheme. We have used two different sets of microscopic interactions: the Argonne v' s implemented 
with the Urbana IX three-body force, and the truncated Urbana V14 plus the Urbana VII three-body force. Binding 
energies have been presented in Tab. IIIII and this is the main result of our work. With the noticeable exception of 
the 12 C nucleus, our binding energies are about 3-4 MeV/nucleon above the experimental ones, consistently with our 
previous results and with those obtained in nuclear matter. 

In our calculations we did not include the contributions of the spin-orbit terms of the correlation and of the p> 
8 terms of the two-body potential, which in 0] have been estimated by means of a local density interpolation of 
the nuclear matter results. These contributions are attractive and of the order of 1 MeV/nucleon. This produces 
differences with the experimental values of the order of those presented in where only the 16 O and 40 Ca nuclei 
were considered. A number of possible sources of the remaining disagreement should be investigated such as the 
three-body correlations and perturbative corrections of the two-body correlations. 

The 12 C nucleus, strongly underbound, stands out of this general trend. As we have already said some specific 
characteristics of this nucleus, not relevant in the other nuclei we have considered, is not properly treated in our 
calculations. We think this could be the presence of deformation phenomena. 

Short-range correlations produce small effects on the one-body density distributions. The distributions at the center 
of the nucleus are lowered by the correlations. This effect is relatively large when only scalar correlations are used, 
and it is strongly reduced when the other terms of the correlations are included. These results are in agreement with 
the findings of 29] where a simple first-order expansion of the density has been used. 



10 



The study of the two-body densities show that the short-range correlations strongly reduce the probability of 
finding two nuclcons at distances smaller than 0.5-0.6 fm. The scalar and isospin, two-body densities, extend on the 
full nuclear volume, while the other densities have a maximum range of about 3 fm. 

Our work shows that microscopic variational calculations for medium-heavy nuclei have reached the same degree 
of accuracy as nuclear matter calculations. Our approach provides a good description of the short-range correlations, 
but it shows some deficiency related to the proper inclusion of long-range correlations. The next obvious step would 
be to extend the theory to have a better treatment of these last type of correlations. 
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l '2 n 


1 fi 

IbQ 


4U Ca 


48 Ca 


2U8p b 


v's + 


d c 


1.20 


2.10 


2.15 


2.10 


2.20 


UIX 


d t 


3.30 


3.70 


3.66 


3.70 


3.60 


via + 


d c 


1.40 


2.10 


2.15 


2.10 


2.20 


UVII 


dt 


3.30 


3.80 


3.86 


3.90 


3.80 



TABLE I: Values, in fm, of the healing distances, obtained minimizing the energy functional J2J with the interactions we have 
used in our calculations. 









ib 


4U Ca 


48 Ca 








1.000 


1.000 


1.000 


1.000 


0.999 


s? 


(SOC) 


0.997 


1.006 


1.008 


1.013 


1.024 


ST 


(J) 


1.000 


1.000 


1.000 


0.999 


0.999 


ST 


(SOC) 


0.997 


1.006 


1.008 


1.021 


1.027 


s 2 


(J) 


1.004 


1.003 


1.001 


1.000 


0.998 


s 2 


(SOC) 


0.996 


1.025 


1.022 


1.042 


1.055 



TABLE II: Normalization of the one- and two-body densities Si and 5*2 respectively. For the one-body densities we give 
separately the contributions of the protons (p) and neutrons (n). The label J (Jastrow) indicates the calculation done with the 
scalar, fx, correlation only, while the SOC labels indicate the results of the complete calculation. 









ib 


4U Ca 


48 Ca 






T 


27.13 


32.33 


41.06 


39.64 


39.56 




v 2-body 


-29.13 


-38.15 


-48.97 


-46.60 


-48.43 




V LS 


-0.25 


-0.38 


-0.39 


-0.35 


-0.45 


v' s + 


Vcoul 


0.67 


0.86 


1.96 


1.57 


3.97 


UIX 


T + V(2) 


-1.58 


-5.34 


-6.34 


-5.74 


-5.35 




Vi-body 


0.66 


0.86 


1.76 


1.61 


1.91 




E 


-0.91 


-4.49 


-4.58 


-4.14 


-3.43 




T 


24.63 


29.25 


37.70 


36.47 


36.48 




V 6 

v 2-body 


-27.08 


-35.84 


-47.16 


-44.86 


-46.87 




Vls 


0.05 


0.03 


0.07 


0.09 


0.04 


via + 


Vcoul 


0.68 


0.88 


2.02 


1.59 


4.03 


UVII 


T + V(2) 


-1.72 


-5.68 


-7.37 


-6.71 


-6.32 




Vi-body 


0.54 


0.69 


1.28 


1.15 


1.41 




E(UIX) 


-1.18 


-4.99 


-6.09 


-5.56 


-4.91 




Eexp 


-7.68 


-7.97 


-8.55 


-8.66 


-7.86 



TABLE III: Energies per nucleon in MeV. We have indicated with T the kinetic energy, with V^-body t ne contribution of the 
first six channels of the two-body interaction, with Vls the spin-orbit contribution, with Vcoul the contribution of the Coulomb 
interaction and with V^-body the total contribution of the three-body force. The raws labeled T + V(2) show the energies 
obtained by considering the two-body interactions only. The experimental values are taken form |23]. 



12 




FIG. 1: Two-body correlation functions f p , as a function of the internucleonic distance. The symbols indicate the correlation 
functions obtained for symmetric nuclear matter (see text for explanations). 
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FIG. 2: Proton density distributions. The full lines are the mean field distributions, the dotted ones have been obtained by 
using scalar correlations (/i) only, and the dashed lines show the results of the complete calculation. 




FIG. 3: The same as Fig. [5] for the neutron density distributions. 
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FIG. 4: Normalized ratios [12( iti for proton and neutron distributions. The meaning of the lines is given in the legend, where we 
have indicated with J the results obtained with scalar correlation (/i) only, and with SOC the results of the full calculation. 
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FIG. 5: Normalized ratios 1271 for charge density distributions. The meaning of the lines is given in the legend. 



15 




FIG. 6: Proton-proton scalar (p=l) two-body distributions l|29fl as a function of the internucleonic distance. The dotted lines 
show the uncorrelated two-body distributions calculated by considering only the first term of eq. 1121 , the direct one. The full 
lines show the total uncorrelated two-body distributions as given by eq. Ill 2t . The other two types of lines indicate the results 
obtained with scalar correlations only (Jastrow) and in the complete calculation (SOC). 




FIG. 7: The same as Fig. E]for the proton- neutron case. The dotted and full lines coincide, since the exchange term of eq. 
1121 1 is zero. 
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FIG. 8: The same as Fig. E]for the neutron-neutron. 




FIG. 9: Proton-proton and neutron-neutron two-body distribution functions 12911 of 208 Pb nucleus. As discussed in the text, 
for equal particles, the p=2,4,6 distribution functions are respectively equal to the p=l,3,5 ones. 




FIG. 10: The same as Fig. E]for the proton-neutron case. 



